Meridional overturning circulation: stability and ocean feedbacks in a box model. 
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A box model of the inter-hemispheric Atlantic meridional overturning circulation is developed, 
including a variable pycnocline depth for the tropical and subtropical regions. The circulation is 
forced by winds over a periodic channel in the south and by freshwater forcing at the surface. The 
model is aimed at investigating the ocean feedbacks related to perturbations in freshwater forcing 
from the atmosphere, and to changes in freshwater transport in the ocean. These feedbacks are 
closely connected with the stability properties of the meridional overturning circulation, in particular 
in response to freshwater perturbations. 

A separate box is used for representing the region north of the Antarctic circumpolar current 
in the Atlantic sector. The density difference between this region and the north of the basin is 
then used for scaling the downwelling in the north. These choices are essential for reproducing 
the sensitivity of the meridional overturning circulation observed in GCMs, and therefore suggest 
that the southernmost part of the Atlantic Ocean north of the Drake Passage is of fundamental 
>— ^ importance for the stability of the meridional overturning circulation. With this configuration, the 

magnitude of the freshwater transport by the southern subtropical gyre strongly affects the response 
^ of the meridional overturning circulation to external forcing. 

The role of the freshwater transport by the overturning circulation (M ov ) as a stability indicator 
is discussed. It is investigated under which conditions its sign at the latitude of the southern tip 
of Africa can provide information on the existence of a second, permanently shut down, state of 
the overturning circulation in the box model. M ov will be an adequate indicator of the existence 
of multiple equilibria only if salt-advection feedback dominates over other processes in determining 
the response of the circulation to freshwater anomalies. M ov is a perfect indicator if feedbacks other 
than salt-advection are negligible. 
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I. INTRODUCTION 

The stability of the Atlantic Meridional Overturning Circulation (MOC) in the ocean has been a topic of invest- 
igation in oceanography and climatology since the pioneering work of Stommel pQ. The central idea of Stommel 
PQ is that the MOC is linked to a feedback between the circulation and the salinity advected from the tropics to 
the subpolar regions, increasing density where sinking occurs and enhancing the MOC strength. This mechanism, 
usually referred to as salt-advection feedback, is also responsible for the collapse and reversal of the circulation, if 
a sufficiently strong freshwater anomaly is dumped in the subpolar areas. This concept was originally developed to 
describe the MOC in a single hemisphere, with the circulation being driven only by buoyancy forcing and, implicitly, 
by mixing-induced vertical diffusivity at low latitudes. It has been used to describe the inter-hemispheric MOC 
as well, with the circulation being driven by a north-south, instead of equator-pole density difference [2]-[4]. Even 
if extremely simple, this paradigm has been helpful in understanding the behaviour of the MOC in more complex 
General Circulation Models (GCMs) of various complexity, which all show that the MOC can substantially weaken if 
a sufficiently strong perturbation is applied. 

It has been shown in several earth system models of intermediate complexity that if the freshwater forcing over the 
North Atlantic (say measured by a parameter 7) is increased, the MOC will first slow-down and then undergo a more 
or less abrupt transition to a collapsed state [5] . The point of collapse is believed to be associated, in the terminology 
of dynamical systems theory, to a saddle node bifurcation. At such a point, which we will call L\, a previously stable 
steady state of the system (the present day "ON" state of the MOC, see Fig. [I]) loses its stability, and the system 
jumps to another stable steady state with little or reversed overturning in the Atlantic (the "OFF" state, see Fig. [I]). 
If the forcing is then reduced again, the MOC does not recover immediately, but stays collapsed until another critical 
value of the freshwater forcing is reached. At this second point, named I/2, the "OFF" state loses its stability and 
the system jumps to the "ON" state, with vigorous overturning. In the region of 7 between L\ and L2, called the 
multiple equilibria (ME) regime, two stable steady states coexist under the same boundary conditions, separated by 
an unstable steady state of the system, usually not observable in numerical models performing time integration. 

In this context, it has been suggested that the freshwater transport by the MOC at the southern border of the 
Atlantic Ocean, usually shorthanded s , may play a special role in determining MOC stability [H[6]. In analogy 
with the original salt-advection feedback of the Stommel model, the sign of the freshwater transport may determine 
whether the MOC could be permanently collapsed or, in other words, if a second stable steady solution, with a very 
weak or reversed MOC, exists. 

If s is negative, the MOC is importing salt into the basin. A negative perturbation on the MOC would 
determine a negative salinity anomaly, and consequently a negative density anomaly, within the Atlantic basin. If 
sufficiently strong, this density anomaly can be advected northward up to the sinking regions of the northern North 
Atlantic, either by the gyres or by the MOC itself. There, it would feed back onto the MOC strength, reducing 
the downwelling rate and amplifying the initial negative MOC perturbation, providing a mechanism for the collapse 
of the MOC. If, on the other hand, s is positive, the MOC is importing freshwater into the Atlantic basin. A 
negative perturbation in MOC strength would then lead to a positive density perturbation, as it would bring a positive 
salinity anomaly in the Atlantic, and no amplification of the initial perturbation is possible. In this second case, even 
if the MOC can be slowed down substantially or even reversed by a sufficiently strong freshwater perturbation, it 
will spontaneously recover when the perturbation is removed, as no stable "OFF" state of the MOC exists under the 
unperturbed boundary conditions. 

In de Vries and Weber [Hj , it is shown that bringing s to negative values through appropriate surface freshwater 
flux anomalies it was possible to permanently collapse the MOC if net evaporation over the North Atlantic is below 
a critical value. In particular, it was shown that a dipolc freshwater anomaly applied over the southern portion of 
the Atlantic Ocean affects the freshwater transport by the southern subtropical gyre circulation (usually referred to 
as M^°° s ), and hence can affect Af30°S. with an increase of ^? ^ can eventually become negative, enabling 
a permanent shut-down of the MOC. 

Recently, a series of numerical experiments with two different GCMs was performed to further explore the import- 
ance of freshwater budget of the Atlantic Ocean. In Cimatoribus et al. 7 , it is shown that the freshwater transport 
by the southern subtropical gyre is of paramount importance in determining the stability of the MOC, even if the 
overturning rate is only weakly sensitive to s . Similarly to de Vries and Weber [5J, s and s were 
tuned through freshwater anomalies at the surface, showing that the sensitivity of the MOC to changes in the net 
precipitation over the Atlantic Ocean depends strongly on how the freshwater transport at the southern border of the 
Atlantic is divided between the gyre and the MOC. In particular, as the freshwater transport by the gyre is increased, 
Mj^j s decreases and this leads to a permanent collapse of the MOC for weaker reductions of the net evaporation. 
The MOC can even be collapsed by solely increasing M^° s , without changing the integrated net evaporation over 
the Atlantic basin. A similar sensitivity of the MOC to the freshwater transport in the South Atlantic can be inferred 
from the work of Marsh et al. [S]. They show that, in a low resolution model, changes in the zonal salinity contrast 
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over the South Atlantic control the MOC sensitivity to freshwater perturbations. 

The results above consistently point to the importance of the freshwater (that is, buoyancy) forcing in setting the 
MOC stability and strength in connection with the salt-advection feedback, but they are in apparent contradiction 
with the point of view of Ferrari and Wunsch [lOj . who suggest that the direct mechanical forcing by the wind over the 
Southern Ocean is the main driver of the MOC. A synthesis of these two views was attempted in Wolfe and Cessi [TT] 
and Nikurashin and Vallis [12]; a series of numerical experiments with an ocean-only GCM (both at eddy-resolving 
and non eddy-resolving resolutions) contribute to a description of the MOC as driven by the wind-induced upwelling 
in the Southern Ocean, but with buoyancy distribution controlling whether an active pole-to-pole overturning can 
actually take place. A similar synthesis of these two views was attempted in simple models by Gnanadesikan |14j and 
Johnson et al. [13] (from now on J07). 

On the other hand, the results of Longworth et al. [9] point to the importance of the transport of salinity by the 
wind driven gyres (represented in a box model by horizontal diffusion) for the stability of the MOC. In particular, 
they show that a simultaneous increase in the horizontal diffusivity in both hemispheres is in general a stabilising 
factor for the MOC, providing a way to transport salt from low to high latitudes that bypasses the salt-advection 
feedback mechanism. 

In this paper, the main goal is to understand the results of Cimatoribus et al. [7J in the framework of a box model, 
and to better understand the validity of s as an indicator of MOC stability. The box model is based on that 
developed by J07 and Gnanadesikan [14) . with a few important differences. Its formulation and solution method 
are described in Sec. [TTJ We are particularly interested in understanding the role of the freshwater transport by the 
gyre and overturning circulation at the southern entrance of the Atlantic Ocean on the changes in stability of the 
MOC found in Cimatoribus et al. [7]. This can be done by including a representation of the transport by the southern 
subtropical gyre as a horizontal diffusive transport, similarly to what was done by Longworth et al. !9j. The sensitivity 
of the solutions of the model to this as well as to other parameters of the system is described in Sec. |III| In relation 
to this, the conditions under which s can provide useful information on the stability properties of the MOC are 
discussed in Sec. IIII CI 



II. METHODS 



A. Model definition 



The box model (see Fig. [2]) used in this work consists of a basin spanning both hemispheres (which represents the 
Atlantic Ocean) connected in the south to a periodic channel (representing the Southern Ocean, shorthandcd s). It 
is assumed that the entire flow at depth from the Atlantic Ocean is upwelled into the mixed layer represented by 
box s, due to Ekman pumping. Upwelling in other basins, due to diapycnal mixing, is thus neglected in this work. 
The basin has four other boxes: a pycnocline (made of two boxes: t and ts) and a deep box (box d) separated by a 
pycnocline at a variable depth D, and a northern box (box n) of fixed volume. The northern box only represents that 
portion of the Northern Seas where downwelling actually takes place. Differently from J07 and Gnanadesikan [14] . 
the Atlantic thermocline and the Southern Ocean are connected through a further box, named ts. The box extends 
south of the latitude of the tip of the African continent (approximately 30°S). The use of this box, motivated in detail 
in the following sections, enables to compute the meridional density gradient within the Atlantic basin, assumed 



to be the driver of the overturning circulation (see Sees. II A 1 and II A 2). This region is characterised by larger 



isopycnal slopes than within the Atlantic basin (see, e. g. Lumpkin and Speer |T5|), with the pycnocline becoming 
shallower moving poleward. The box ends in the south at the latitude where the deepest Antarctic Intermediate 
Waters (northward flowing) outcrop, approximately 40° S [T5]. The depth of this box at its northern end is the same 
depth scale D used for the main thermocline box t. This depth scale is obtained from the pycnocline model used 
by J07 and Gnanadesikan [TT], but with the thermocline box divided in the two boxes t and ts, as described. The 
temperature is fixed in all boxes, while the salinities are prognostic variables of the model. This is a common choice 
in simple models of the MOC, stemming from the much shorter decay time of large scale temperature perturbations 
in the ocean, with respect to salinity ones. This choice implies that we will neglect any feedback with the atmosphere. 
This appears as an obvious limitation of the work; but it may not be, as recent studies suggest that atmospheric 
feedbacks play a secondary role for MOC stability [T51U8] . 

The pycnocline model results from a balance between the following terms: 

1. Ekman transport into the Atlantic basin from the Southern Ocean (and through the ts box), associated with 
the wind-driven upwelling in the Southern Ocean (qEk)', 

2. eddy induced flow out of the Atlantic basin (from box ts to box s) due to baroclinic instability (g e ); 



3. diffusively controlled upwelling of water from the deep layer (qu)', 



4. downwelling in the North Atlantic (qn)- 

Within box t, we assume that the MOC can be approximated as a coherent flow, northward in the pycnocline and 
southward at depth. This pycnocline model can be paralleled to the view of the MOC by Wolfe and Cessi [TJJ, which 
also includes a wind-driven upwelling ("pushed" from the south) and a diffusively controlled upwelling inside the 
basin, the latter compensating only for a small part of the downwelling in the north. The volume of the Atlantic 
pycnocline is the sum of the volumes of boxes t and ts, and the volume budget of the pycnocline thus reads: 



In the equation, r is an average zonal wind stress amplitude over the periodic channel region, L x $ is the zonal extent 
of the Southern Ocean, L x a is the zonal extent of the Atlantic Ocean at its southern end, D is the pycnocline depth, 
Po is a reference density, /<j is the Coriolis parameter in the frontal region of the Southern Ocean, Aqm is an eddy 
diffusivity, L y is the meridional extent of the frontal region, A is the horizontal area of the Atlantic pycnocline and 
k is the vertical diffusivity. The total volume flux at the southern boundary of boxes t and ts is called qs; it is the 
difference between the Ekman inflow and the outflow due to baroclinic instability, while qn, the downwelling flux 
in the north, is left unspecified for the moment. As the Ekman inflow in the Atlantic Ocean is related to the wind 
forcing over the whole ACC [THIHO]) the zonal width of the entire Southern Ocean is used in the computation of qEk- 
The zonal width of the Atlantic Ocean is used instead for computing the eddy outflow into the periodic channel. The 
volume of the box t is given by V t = AD, while that of box ts by V ts = L x aL v D /2. It can be seen from Eq. (JTJ) that 
the volume flux from the box s goes through box ts into the main thermocline without any change, we thus neglect 
any diapycnal upwelling within box ts. As far as the pycnocline model is concerned, boxes t and ts behave as a whole, 
and the model is in fact completely equivalent to that of J07 and Gnanadesikan [14] . 

The model can be interpreted as an analogue of the adiabatic pole-to-pole circulation paradigm of Wolfe and Cessi 
[TT] and Nikurashin and Vallis [T2]. Diabatic fluxes mainly take place in the north, at the interface between boxes t 
and n, where thermocline waters are made denser by the interaction with the atmosphere, and in the Southern Ocean, 
where water flowing northward in the Ekman layer gains buoyancy. In the interior, the flow is mainly adiabatic, and 
the lower branch of the overturning is virtually isolated from the upper one, apart from the diffusive upwelling qu, 
very small in our configuration. The upwelling in the south is also adiabatic, as water mass transformation only takes 
place when the lower branch of the MOC enters the box s, the mixed layer of the ACC. 



Based on geostrophy, J07 used the scaling q^ oc D 2 (pd — p t ) /po for the downwelling jT4j [2TJ [22] , with pd the density 
of the deep box and pt that of the thermocline box. Since at steady state the densities of the deep and northern 
boxes are equal, this scaling can be written as qjy oc D 2 (p n — pt)/po, with p n the density of the northern box. The 
density difference between the northern and pycnocline boxes is almost entirely determined, in the box model, by 
the temperature difference (fixed in J07 and here). Using this scaling, the sensitivity to freshwater fluxes is thus too 
low (more than 1 Sv is needed to collapse the MOC — not shown). Furthermore, if a representation of the southern 
subtropical gyre is included in the box model of J07, the MOC shows no sensitivity to the gyre strength and to the 
freshwater transport associated with it (not shown). This latter result is in striking disagreement with the results 
of Cimatoribus et al. [7] , as discussed in section [I] 

Levermann and Fiirst [23J tested the scalings used in J07 and Gnanadesikan [JJ]. They show that the downwelling 
can be described by q?j = rjD 2 Ap, with Ap the density difference between north-western and tropical water at the 
base of the thermocline. Even if the validity of the above results for different models is still uncertain [24], the 
connection between overturning strength and meridional pressure gradients is well established (see e. g. Sijp et al. 
[2"5] . Cessi and Wolfe [35]). In the simple context of the box model, this translates into a meridional density gradient 
in the thermocline |21| . We exploit this idea by taking Ap — (p n — p ts )/po- In our model, the box ts is the southern 
end of the thermocline, which means that the overturning circulation is proportional to the density difference between 
the sinking regions in the north and the southern end of the thermocline. The pycnocline model then reads: 



d(V t + V ts ) = tL xS 
dt p \f s \ 




(1) 



= qEk - qe + qu - qN 

= qs + qu - qN, with q s = qEk - q e - 



1. The scaling for 



5 



From a heuristic point of view, this scaling for g^v implies that, when a large meridional density gradient is present, 
the sensitivity of the MOC is mainly determined by changes in pycnocline depth (neglecting the variable pycnocline 
depth leads in fact to unrealistic sensitivity to wind stress, not shown), while a freshwater flux in the north leads to a 
collapse due to the reduction of the inter-hemispheric density difference. We thus view the collapse due to a freshwater 
perturbation as "buoyancy driven" , even though the sensitivity of the MOC in the "ON" state is dominated by "wind 
driven" dynamics (e. g. changes in the Ekman inflow). Even though we can collapse the MOC exploiting a buoyancy 
flux, reducing the wind stress over the Southern Ocean below a critical value also brings q$ to zero (see Sec. Ill A 1 1, 
leaving a purely diffusive intrahemispheric MOC (qn — qu > 0) similarly to what is suggested by the results of Wolfe 
and Cessi [IT] . 



2. The role of the box ts 



The pycnocline model of Eq. ([2]) is equivalent to the one defined in J07, the only difference being our thermocline 
box split into two boxes, in order to emphasise the importance of the meridional density gradient. 



Considering the salinity transport, instead, we proceed in a different way from J07. We argue that at the southern 
border of the thermocline (the southern border of box ts in our model) the Ekman inflow and eddy outflow, whose 
difference gives the net volume flux, have to be taken into account separately concerning the salt transport. The salt 
flux from box s to box ts has then to be written as qEkS s — q e St s - In other words, even if the net volume flux and 
net density flux may sum to zero, there can still be exchange of salt between the Southern Ocean and the pycnocline. 
In case of no net volume flux (qEk = Qe) the salt transport must be compensated in density by heat transport (not 
considered explicitly in this work). This salinity transport must then be associated to transport of spiciness: it 
represents in this case an isopycnal salinity flux. In general, for q^k 7^ 9e ; the net volume flux will be associated with 
density transport on top of this isopycnal salinity transport. Although the Ekman and eddy fluxes may compensate 
globally, they will not do so locally, so that transport of properties is still present. From box ts to box t, instead, we 
assume that the flow is coherent because the flow is already within the pycnocline layer, in geostrophic balance and 
with less prominent baroclinic eddy flow. The transport of salt from box ts to box t can then be written simply as 
QsSts (for the case of a northward net flow). Upwclling in the box s is also taken proportional to the net flow qs- 



A last connection between boxes ts and t is added through a diffusive constant, to represent in the simplest possible 
way the transport by the southern subtropical gyre, similarly to what was done by Longworth et al. [5j and J07 (the 
latter only considered the northern subpolar gyre). We argue that, if the MOC is sensitive to the difference between 
the sinking regions and the thermocline waters north of the ACC, the southern subtropical gyre may play a very 
important role in setting this density difference. The argument in J07 of the effect of r$ being negligible, at least in 
comparison with the gyre in the northern part of the basin, does not seem to be applicable in this context. 



With this configuration, we neglect the interaction with the atmosphere of the area represented by box ts. This 
choice is done for simplicity, but seems to be at least indirectly supported by Garzoli and Matano |27j , who show that 
little mass water transformation takes place for lighter water classes in this region (their Fig. 6). 



A more important assumption in this work is the neglect of the exchanges between the Atlantic Ocean and the 
other Ocean basins, in particular with the Indian Ocean through Agulhas Leakage. This is still a highly debated 
subject, outside the scope of the present paper. 
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3. Salt conservation equations 

The equations needed to close the system are given by salt and volume conservation in each box. They comprise 
the following set of equations for a state of the MOC similar to the present day one (qN > qs > 1u > 0): 



— -jT — = QsSts + quSd - qNS t + r s {S ts - S t ) + r N (S n - S t ) + 2E 8 S , (3a) 

Vt °f ts ^ = q Ek S s - q e S ts - q s S ts + r s (S t - S ts ), (3b) 
at 

v n ^r = - 5 ») + r *( 5 * - 5 ») - ( E > + E ^ s °> ( 3c ) 

= qNS„ - quS d - qsSd, (3d) 

V s -^ = qsS d + q a S t s - qEkS s - (E s - E a )S Q , (3e) 

SoVtot - V n S n - V d Sd - V t St - V ts S ts - V S S S = 0, (3f) 



with V being the volume of the various boxes and V tot = V n + Vj + V ts + Vd + V S . The last equation states salt 
conservation in the model. S is the salinity in the various boxes, r$ is the gyre exchange between boxes t and ts, 
while r^v is the gyre exchange between the pycnocline and northern box. The salt transport by the gyres takes in 
this simple framework the form of a diffusive constant. E s is the symmetric part of the atmospheric freshwater flux, 
from the pycnocline to the boxes n and s. E a is the asymmetric part of the atmospheric freshwater flux, transferring 
freshwater from the Southern Ocean into the box n. Physically, the asymmetry in the freshwater forcing is due to 
the different sizes of the boxes in the north and in the south, differences in precipitation amount connected with the 
fact that box n is closer to the pole than box s and differences in the runoff from the continents. We will change 
E a , changing the North-South density gradient, for collapsing the MOC. In reality, possible sources of the freshwater 
flux are ice-sheets at high latitudes or Arctic sea-ice, not represented in the box model; these will generally cause a 
net change in the total freshwater content of the ocean. Using E a as a "hosing" flux, enables to maintain a closed 
salinity budget, and thus to explore steady states of the system; the approach is similar to that used in GCMs, where 
"hosing" fluxes are usually compensated over the rest of the Ocean surface. Both E s and E a are represented through 
a virtual salt flux. Despite the different definition used here, the atmospheric freshwater fluxes E a and E s can be cast 
in a form equivalent to that used in Longworth et al. [9J and Scott et al. [18]. So is a reference salinity and a linear 
equation of state is used, i. e. p = Po(l — a(T — T ) + j3(S — So)), where a and /3 are the constant thermal expansion 
and haline contraction coefficients respectively, To is a reference temperature. 

To summarise, there are six differences in our model definition from J07: (i) a different choice of the box sizes, 
to account for our interpretation of the model as the Atlantic and Southern Oceans alone; (ii) the pycnocline layer 
is split in two meridionally at the latitude of the southern tip of Africa; (iii) the two components of salt advection 
between the southern box and the thermocline, due to the Ekman and eddy fluxes, are treated separately; (iv) the 
meridional density gradient used in the scaling of qjv; (v) the use of a term representing gyre exchange (rs) between 
boxes t and ts and (vi) a slightly different parameter choice including lower vertical and eddy diffusivities. 



4- Reference solution 

At steady state, if q^ > qs > 0, it follows that Sd = S n from Eq. (|3d[ ), as qs = qN — qu- The reference parameter 
values are shown in Tab.|T] A solution for the system of Eqs. Q and (|3J) can be found using Mathematica software [28], 
but only after substituting the parameters with their numerical values. The solution for the reference configuration 
of Tab. |T]is reported in Tab. |TT} 

Similarly to J07, the reference solution consists of a flow northward in the boxes ts and t and southward at depth, 
downwelling from box n to box d and upwelling from box d to box s. Diffusive upwelling takes place from box d 
to box t. The northern box is saltier than the southern one under symmetric forcing, due to the asymmetry of the 
circulation itself (in other words, the circulation is maintained by the salt-advection feedback). Differently from J07, 
the value of pycnocline depth is close to what commonly observed in the real ocean and in GCMs. This is mainly due 
to smaller downwelling flux than in J07, closer to the present day estimate of 18.7 ± 2.1 Sv [29] . given by the different 
scaling for q^ and by the smaller vertical diffusion considered here. The value of the hydraulic constant 77 used is 
the one estimated in [23], but we stress that this quantity is badly constrained since it will strongly depend on the 
choice of the temperature difference between boxes n and ts (or, in a more realistic framework on the locations used 
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for the computation of the meridional density difference) . For this reason, we will study the model at different values 
of the hydraulic constant r\ (0.5, 1 and f .5 times the value given by Levermann and Fiirst [23]). in order to assess the 
robustness of the results. 

The temperature chosen for the box T ts is 10°C, giving a difference of 5°C between the northern box and the 
southern end of the pycnocline. This temperature is representative of the waters north of the periodic channel in 
the South Atlantic. This temperature, or the temperature difference between boxes n and ts |37] . is particularly 
important since it determines the strength of the MOC for fixed rj and a given salinity distribution. The choice of 
this temperature difference was made in order to use the estimate of rj given by Levermann and Fiirst j23j , the most 
recent we are aware of, and to obtain at the same time a realistic MOC rate. 

A negative value of E a is used for the reference solution, to account for the different interpretation of the boxes in 
the model, with respect to J07. A negative E a is needed in order to reduce the export of freshwater from the Atlantic 
Ocean, and also to reduce the net precipitation over the sinking regions of the North Atlantic, also considering the 
asymmetry of the sizes of box n and s. 



5. Freshwater budget 

In order to compare our results with those from more complicated models, we compute the equivalent freshwater 
budget of the "Atlantic" basin north of 30°S (see Fig. [2]). With the definitions: 

M ov = --^qs (S ts - S d ) , 

DO 

M az = -^r s (S ts -S t ), 

DO 

valid under the same conditions considered in Eqs. ([3]), we find that the net steady evaporation from the Atlantic 
basin must be balanced by the freshwater import by the ocean circulation: 

(E s - E a ) = M ov + M az . (4) 

This budget includes the dominant terms in GCMs of various complexity [71|3D]. Eq. Q states that the freshwater 
transported through the atmosphere out of the basin must be balanced by the freshwater transport by the ocean 
circulation, split into its overturning (M ov ) and gyre (M az ) components. M az can be controlled by varying rg, and 
M ov will change in the opposite way, as long as the net evaporation is kept constant. This behaviour, which is trivial 
in this model, has been observed in GCMs as well. In Cimatoribus et al. [TJ, M az is changed using anomalies in the 
surface freshwater flux instead of perturbing the gyre strength, and M ov reacts in a similar way as here, compensating 
for the changes in the azonal transport and keeping the total budget closed. 



B. Solution method 



The system of equations ^ and Q is studied using the software AUTO-07p, which enables the continuation and 
bifurcation of solutions of systems of Ordinary Differential Equations (ODEs) 31J. Only exact (within numerical 
accuracy) steady states are obtained with this technique, and no transient behaviour is studied in any part of this 
work. 

The system of equations ^ and ^ can be written as the autonomous ODEs system: 

^ = G(x W ,p), (5) 

with x the 6-dimensional state vector, p the vector containing the system parameters, and G a nonlinear mapping. 
Given a steady state solution xo of ^ for a particular parameter set po, AUTO-07p can track the evolution of the 
steady-state solution as the value of one parameter of the system is changed; the solutions can be stable or unstable, 
and their stability is determined as well. The initial solution x can be analytical or numerical, as long as it is known 
with sufficient numerical precision. The software is based on Newton's method for continuation of steady states; for 
further details, see Doedel and Oldeman [3T] and references therein. Special points, such as saddle node bifurcations, 
are detected and can be continued in two parameters, obtaining the locus of those special points in a parameter 
plane. Multiple steady states can be identified following the solution beyond a special point and, more in general, a 
systematic exploration of the parameter space is possible. All the results are obtained with small continuation step 
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sizes, so that the results are not sensitive to the value actually used. The relative convergence criterion for the steady 
solution in AUTO-07p is set to KT 7 . 



In the next section, starting from the reference solution obtained from Mathematica (see Sec. II A 4), the sensitivity 
of the model to several parameters is studied. Solutions are continued within physically plausible regions of the 
parameter space, or down to the point where qs = 0. When the sensitivity to freshwater fluxes is studied, the solution 
is continued for qs < down to the point at which the downwclling in the north stops. When qs changes sign, the 
definition of the box model has to be changed. In particular, Eqs. (3a), (3b), (3d) and (3e) become for a negative qs- 



d(VtS t ) 
dt 

d {VuSt.) 
dt 

d(V d S d ) 



dt 



dt 



qsS t + quSd - qNS t + r s (S ts - S t ) + r N (S n - S t ) + 2E S S , 

qEkSs - q e S t s ~ qsS t + r s (S t - S ts ), 

QnSu — qi/Sd — qs$s, 

qs$s + q e S t s — QEkS s — (E s — E a )S , 



while the equations for the northern box and total salt conservation are unchanged. The two sets of equations are 
automatically chosen in order to match the solution at qs = 0. 

The solutions are not continued below q^ — as then the MOC is completely reversed and the scaling for q^ 
is not meaningful anymore. In fact, if qw was allowed to reach negative values, the scaling used in the "ON" state 
would represent an enhanced upwelling in the high latitudes of the North Atlantic, for which no plausible physical 
mechanism is known. Furthermore, a reversed MOC is determined by the shallow outflow from the Atlantic, in the 
present box model this is due to an eddy flux stronger than the Ekman transport. Whether such a situation could 
exist in reality is unclear. With these limitations in mind, the study stops at the point where q_/v = with qs small 
and negative. The change in sign of qs when the downwelling in the north is weaker than the mixing-driven upwelling 
((/TV < Qu) does not affect the fundamental properties of the model, which is still controlled by the same scaling law for 
gjv- On the contrary, the model is no more valid when Ap and q^ change sign. An extension of the model including 
a consistent representation of the completely reversed MOC is beyond the scope of this work. 



III. RESULTS 



This section is comprised of three subsections. In subsection |III A| we investigate the sensitivity of the model 
solutions to the southern ocean wind stress and the strength of the southern gyre. In subsection III B the stability of 
the MOC to freshwater perturbations is considered while in subsection |III C[ the results in |III B| are interpreted using 
the freshwater budget of the box model. 



A. Model sensitivity 

In general, the sensitivity to the model parameters is very similar to that found in J07. We will briefly discuss the 
importance of wind stress, and then analyse the sensitivity of the model to the gyre transport and to changes in the 
freshwater fluxes. 



1. Southern Ocean wind stress 



In Fig. [3] the sensitivity of D, qs and qjy to the wind stress r is shown. When the wind stress is changed, the 
role of pycnocline depth in determining the MOC strength is dominant. The values of the volume transports change 
mostly due to the increased Ekman inflow at the southern border, generating a deeper pycnocline and a consequently 
stronger outflow in the north, which is qualitatively the same behaviour as in J07. The changes in the meridional 
density difference, with a reduction of the salinity in the box ts (due to the increased inflow from the fresher box 
s) but little change in the northern box, enhance the changes in q^ and qs- Below a certain critical value of t 
(approximately 0.02 Nm -2 ) the eddy flux becomes as large as the Ekman inflow, and the waters sinking in the north 
are upwelled only through diffusion within the Atlantic basin. Without a net inflow from the Southern Ocean, the 
MOC can not extend beyond the Atlantic basin, consistently with what is shown in the numerical experiments of 
no pole-to-pole MOC is possible without wind-stress over the Southern Ocean. 



9 



2. Southern gyre 



As discussed in Sec. II A 4 a basic representation of the southern subtropical gyre is included in this work. The 
sensitivity of D, qs and qn to r$ is shown in Fig. [4] In addition, the sensitivity of M ov and M az is also shown in 
the bottom panel. Both qs and q^ decrease when the gyre strength increases, with stronger sensitivity for lower 
values of rg. This is a purely buoyancy driven response, which we could not reproduce using the J07 model with a 
term equivalent to rs added to their equations (not shown). Downwelling in the north decreases due to the lowering 
of the north-south density difference, caused by the increased exchange between boxes ts and t which leads to a 
salinification of the southern end of the thermocline and a slight freshening of the box t. The pycnocline depth 
increases moderately, in connection with the slightly stronger decrease of q^ compared to qs, and a slightly decreased 
upwelling flux qjj. The northward volume transport into box ts is reduced by the stronger eddy outflow, associated 
with a deeper pycnocline. This response is absent from J07 since r$ strongly increases the salinity of box ts, but 
very little that of boxes t, n and d. The scaling used in this work for qpj is thus affecting the MOC (as it depends on 
(Pn ~ Pts)), contrary to the one used in J07 (oc (pj — pt))- 

As expected from the freshwater budget, Eq. Q, M ov decreases and becomes negative as M az increases from zero, 
similarly to what is observed in numerical experiments with GCMs [HIS1[7]. This behaviour, trivial in the box model, 
has been exploited in GCMs to change M ov , in order to reach the ME regime [SI [7]. In GCMs, this behaviour is 
usually induced by different means; M az is increased not by modifying the gyre strength, but rather by perturbing 
the surface salinity, and M ov compensates M az due to changes in the intermediate depth ocean stratification [7]. 

The response of M az to rs saturates for large values of rs- This is due to the decrease in the salinity difference 
between boxes ts and t as rs increases, which limits the ability to tune M az using rs in the box model. Another 
limit of representing the gyre through a diffusive constant, already discussed by Longworth et al. ]9], is the fact that 
the freshwater exchange depends on the meridional density difference while in reality the freshwater transport by the 
wind-driven gyre will rather scale with the zonal salinity gradient. 



B. Stability to freshwater perturbations 

1 . Symmetric freshwater flux 

Opposite to what found in J07, the MOC strength increases when E s is increased (Fig. [5j the results are shown for 
three different values of 77). The different sign of the sensitivity to E s is due to the new definition of the qN scaling. 
Increasing E s , freshwater is moved from the thermocline box to the two high-latitude boxes. Given the sense of the 
overturning circulation, this produces an increase in the density difference between boxes n and ts, as the freshening 
by the freshwater forcing in the box n is partly compensated by increased salt advection from the thermocline, while 
box ts is made fresher by its exchange with the southern box. The response is thus due to the asymmetry between 
north and south, induced by the MOC itself; the sense of the circulation (northward above the pycnocline) causes a 
stronger salt transport from the thermocline towards the north than towards the south. It is interesting to note that 
an increase of E s amounts to a decrease of the freshwater forcing over the entire Atlantic basin, as the net evaporation 
from the Atlantic (boxes ts, t and n) is given (E s — E a ). The pycnocline depth decreases in response to the stronger 
outflow from box t to box n. 

Extending the study to the non physical regime E s < (net evaporation in the high latitudes) confirms the role 
of the MOC asymmetry; the MOC can not be collapsed decreasing E s below some threshold because, as the MOC 
weakens, also the impact of E s on the north-south density difference decreases. The sensitivity of the MOC thus goes 
to zero for large enough negative values of E s (not shown). Similarly, as the value of rs increases, the sensitivity of qN 
to Es decreases (not shown) as the southern subtropical gyre transports salty waters from box t to box ts regardless 
of the sense of the overturning circulation. 



2. Asymmetric freshwater flux 

In Fig. [6j three bifurcation diagrams computed for different values of rj are shown, including D, gjy, qs and M ov . 

As expected, we find that an increase in the freshwater flux E a leads to a decrease of the MOC strength. We thus 
recover the result of [TB] , that the MOC strength is controlled by the net evaporation of the overall Atlantic basin. 
Even if the area where a freshwater anomaly is applied determines quantitative changes in the sensitivity of the MOC, 
the sign of this sensitivity is the same as long as the perturbation is applied within the basin. In the framework of 
our box model, both increasing E a or decreasing E s amounts to a reduction of net evaporation out of the Atlantic 
basin, and both determine a decrease of the MOC strength. An increase of E a is, however, much more effective in 
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reducing the MOC strength than a decrease of E s , since it changes directly the north-south density difference to 
which qjv is proportional, and can completely collapse the MOC, differently from E s , as discussed in the previous 
section. The freshwater flux needed to collapse the MOC is about 0.4 Sv, a value comparable to that used in several 
"hosing" experiments with models of various complexity (e. g. [32 ) and much smaller than the symmetric freshwater 
flux needed to collapse the J07 model. Even if a similar freshwater flux E a is used as a "hosing" in the model of J07, 
the sensitivity in this model is still relatively weak. 

In Fig. [6j the solution is continued to negative values of qs, down to the point where the downwelling stops in the 



north. This implies that, at qs = 0, Eqs. (3a|, (3b), (3d) and (3e 



are changed in order to account for the reversal of 
The continuation to negative qs, down 



IIB). 



the advection direction at the boundaries of the boxes s and ts (see Sec 

to the point qat = 0, is performed for discussing the relevance of M ov as a stability indicator (Sec. IIIC) 

At the point where q^ reaches zero, an increase of the freshwater flux E a maintains the collapsed state with q^ = 
and a weakly negative qs compensating the upwelling in the Atlantic, as an increase of E a decreases p n with respect 
to pts- On the other hand, if starting from the point where qjy reaches zero the freshwater flux is decreased further, 
the reversed solution can no longer exist, as a decrease of E a determines an increase of p n with respect to p ts , which 
must be associated at the steady state with downwelling in the north. We thus conclude that the point defined by 
qN(E a ) = is the second saddle node bifurcation L2 of the bifurcation diagram of the MOC strength of the box 
model (Fig.[l]). 



3. Sensitivity of MOC stability to changes in the gyre circulation 



We now turn to the assessment of the sensitivity of the collapse position, the saddle node bifurcation L\, to changes 
in the gyre circulation. The position of the saddle node bifurcations of Fig. [6] is tracked while rj or rjv are changed. 
This enables to compute the critical value of E a in a range of values of rs (Fig. [7]) and rjy (Fig. [9]), for each value of rj 
used in the bifurcation diagram of Fig. [6] The regime diagram obtained divides the parameter plane in two regions: 
below the critical curve an "ON" state of the MOC is possible, while above the line only a collapsed state is present. 

Considering first the dependence of Li position on rs (top panel of Fig. [7]), we see that the asymmetric freshwater 
flux needed to collapse the MOC is reduced by an increase of the gyre strength in the south. The decrease with rs 
of the freshwater forcing needed for collapsing the MOC is substantial, and leads almost to the disappearance of the 
"ON" state for the highest values of rs and lowest 77. Even if a gyre exchange of 50 Sv or more is unrealistic, the value 
of M az associated with it (Fig. [7j lower panel) is well within the range observed in GCMs [3D] . A more effective way 
of changing M az , i. e. changing the South Atlantic freshwater fluxes, may actually lead to the total disappearance of 
the "ON" state for realistic E a values [7 . Longworth et al. [9] , on the other hand, did not observe this behaviour in 
their box model, finding that the effect of the gyres was always a stabilising one. The difference in results is due to 
the fact that they only considered the impact of changes in both northern and southern gyres at the same time, in a 
much more symmetrical configuration (no periodic channel and ACC in the south, boxes of equal sizes in the north 
and south, and no box ts). 

The sensitivity of the position of L\ can be understood considering the lower panel of Fig. [7J where the value of E a 
at L\ is plotted versus M az instead of versus rg. The loss of stability of the "ON" state is connected to an increase 
of the amount of freshwater transported by the southern gyre into the Atlantic basin. This reduces the north-south 
density difference by increasing the salinity of the box ts and slightly decreasing that of the box t; a smaller E a is 
then sufficient to bring (p n — p ts ) to zero, stopping the sinking in the north. We conclude that an increase of M az , in 
this case obtained by increasing rs, brings an effective negative density perturbation in the basin. In terms of changes 
in the bifurcation diagram, we have a shift to the left of the bifurcation diagram and a slight decrease of the width of 
the ME region as rs increases (see Fig.[8j top panel). 

The impact of the northern gyre on the stability of the MOC (Fig. |9| is opposite to that of the southern gyre. 
A strengthening of the northern gyre increases the salinity in the northern box, increasing the MOC strength (not 
shown) and reducing the sensitivity of the MOC to freshwater perturbations. This translates into two effects: (i) L\ 
moves to higher values of E a , i. e. a stronger freshwater flux from outside the Atlantic into the n box is necessary 
to collapse the MOC and (ii) the collapse of the MOC takes place at weaker MOC rates (Fig. [9J). A strong gyre in 
the north brings salty waters in the n box independently of the state of the MOC, and thus reduces on one hand the 
effectiveness of a freshwater anomaly (£1 moves to higher E a ), and on the other hand the importance of the salt- 
advection feedback {L\ takes place at lower qs and gjv), as it reduces the relative importance of the salt transport by 
the MOC. In connection with the latter point, also the width of the ME regime shrinks markedly (see Fig. [8j lower 
panel). This response to changes in r^ is consistent with what was found by Longworth et al. [3]. 
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C. The role of M ov 

In Fig. [6j the bifurcation diagram of M ov as a function of E a was shown, together with the MOC strength. M, 
monotonically decreases as the freshwater flux in the north is increased until Li is reached. A decrease in M ov is thus 
a robust indication of a reducing distance from the limit point L-y, where the "ON" state is no longer stable. It must 
be kept in mind that all the results presented here refer to steady state solutions. Changes in M ov are thus differences 
between different steady states, and not trends in time. The validity of the box-model approach, and that of M OVl 
far from the steady state is unclear; for this reason we focus on the equilibrium solutions alone. 

We can then consider the ME region, that is the part of the bifurcation diagram between L\ and the end of the 
"OFF" branch, L2 . In this region, a shut down state of the MOC is present besides the "ON" state under the same 
boundary conditions. As discussed in the introduction, it has been suggested that the change in sign of M ov may 
mark the entrance into the ME regime, or in other words M ov may change sign at the value of E a for which L2 is 
reached. If = is taken as L2, we see that M ov is not a perfect indicator of the ME region, but it gives a good 
approximation of the ME region, at least for low values of 77. 

To analyse this issue in more detail, we study the difference between E a \® N , the value of E a for which q^ — and 
Eal^j , the value of E a for which M ov = but qs,qN > 0. The first can be computed considering the system: 

q N = 0, (6a) 

r N {S t - S n ) - {E s + E a )S = 0, (6b) 

qEkS s - q e S ts - qsS t + r s (S t - S ts ) = 0, (6c) 

qsS s + QeSts - qEkS s - (E s - E a )S = 0, (6d) 

E s ~ E a = M ov + M aZ7 (6e) 

qu + qs = 0, (6f) 

where the equations represent, in order, the condition of no downwelling in the north, the salt budget of the northern 
box, the salt budget for the box ts, the salt budget for the box s, the freshwater budget of the Atlantic Ocean (the 
latter three for the case q s < 0, with M ov = q s /So(Sd — S t )), and the volume budget for the pycnocline, all under 
the condition q^ = 0. 

From Eq. (J6f| the pycnocline depth at L2 can be computed; choosing the positive solution of the second order 
equation we obtain: 



^ = ^^(l + Jl + ^^V (7) 



2 L x a Aqm \ \] L y q 2 Ek J 

where we write the value of D at L2 as D\° qN , and the Ekman inflow into the pycnocline, qEk = (tL x s)/ (po \fs\), is left 
implicit. Eliminating the salinities from Eqs. (6a|- ( |6e[ ) and using Eq. @» 

an expression for E a \ can be obtained: 



w ,0 _ P r NqEk~qe\° qN (qe\° qN - qEk + r S ) 

ft I o TV s in / in \ 

r N qsk + q e \° qN {q e \% N - qsk + r s ) 

+ r ATa q e \° N (q e \° qN - q E k + r S ) 

So/3 r N q Ek + q e \° qN (q e \° N - q Ek + r s ) ' 

where q e \q N represents the value of the eddy flux at the point where q^ = and AT = T ts — T n . E a \ qN does not 
depend on 77, as confirmed by the numerical results in Fig. |6] It is interesting to compute the limit of E a \ qN for rg 
and k going to zero. From Eq. Q, we find that: 



lim D\", = 



_ L y q E k 



,qN L t a A 



xA s±GM 



which simply states that if k goes to zero, q E k = <7e if qN — 0, since also qs must be zero. The non diffusive version 
of Eq. Q is then easily obtained: 



k^o yjv r N + r s Sq/3 r N + r s 
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from which the rg — > limit of E a \° can be computed: 



lim E t 



1° 

l?JV 



this shows that the "OFF" solution ends at the point where net evaporation over the Atlantic basin sums up to zero 
(E a = E s ), if no vertical diffusion or gyre transport in the south are present: no "OFF" state is available for values 
of E a lower than E s . It is interesting to note that the presence of r^ is irrelevant in this limit, i. e. the redistribution 
of freshwater within the Atlantic basin is irrelevant. This limiting behaviour is the one discussed in Rahmstorf [3]. 

Similarly, the value of E a associated with M ov = 0, shorthanded E a \° M , can be obtained solving the system: 



Sts — S n 

q N {S t - S n ) + r N (S t - S n ) - (E s + E a )S = 
qEkS s - q e S ts - qsSts + r s (S t - S ts ) = 
qsS„ + q e S ts - qEkS s - (E s - E a )S = 

qs - qN + qu = o 



(10a) 
(10b) 
(10c) 
(lOd) 
(lOe) 



where the first equation is the condition for M ov = if qs > 0, the second is the salt budget for the northern box, 
the third is the salt budget for the ts box, the fourth is the salt budget for the s box (with — S n , as qs > 0) and 
the last equation is the volume budget of the pycnocline. The algebra is in this case more tedious, but an expression 
for E a W lo ^ can be obtained similarly to what done for E a \® N , exploiting Mathematica software [28] (see Appendix |a|. 
The skill of M ov as an indicator of ME can then be measured as 



E„ 



Ea\°M 



The analytical expression obtained for A E and E a \° M are very long, so that they are reported in the Appendix [A| 
The expression of A E in the limit n — > reads: 



lim Ae — — 

k-»0 r N + r s 



ATa 

q Ek ATa -q 



2 V L„ 



1-1/1+4 



AgmL x 



q Ek ATa rj 



~2E S -. j^^^^^^z^^^^^^ 

( te+ r N+ r s )AT an l(^) 2 l-/l + 4(i)' 



q E k ATa rj 



(11) 



from which it trivially follows that —> if also rs —> 0. We thus obtain the important result: 

lim A E = 0, 

rs, k— >0 



(12) 



recovering the result of the simpler model of Rahmstorf [3] . Indeed, M ov is a perfect indicator of ME if the response of 
diffusive upwelling and southern gyre transport to MOC changes are neglected. From Eq. (12) we can expect M ov to 



be a good indicator of the ME regime as long as feedback mechanisms other than the salt-advection feedback, which 
is measured by M ov , do not play a relevant role. In particular, vertical diffusion and the freshwater transport by 
the gyre circulation (or by diffusion) at the southern border of the pycnocline must not provide an effective feedback 
mechanism to changes in the MOC, which in the simple framework of the box model is the limit K,rs — * 0. In the 
limit of no vertical diffusion and no gyre in the south, other mechanisms are irrelevant to the skill of M ov as an 
indicator of ME. 



In Fig. 10 



the dependency of A E on rs, K and rj are shown (from top to bottom). The top panel shows the 
dependency of A^ on rs in the non-diffusive limit n — > and the centre panel that of A^ on k in the limit of 
rs — > 0. In the case of the southern gyre, A^ first grows towards negative values, reaching a minimum at about 5 Sv 
to increase again and reach positive values for higher rs values. An increase in rs has two effects: (i) it increases M az 
and consequently decreases M ov (see Fig. |ZJ bottom), making E a \° M smaller, and (ii) it shifts Li to lower values of 
E a (see Fig. [8j top). The first change depends on how effective rs is in increasing M az , and from Fig. |i]it is clear that 
M az increases quickly for low rs, but then saturates for larger values of rs, as the salinity difference between boxes t 
and ts decreases. The value of E a \° qN decreases faster than E a \\ { at small rs, but then saturates as well being almost 
constant above rs ~ 30 Sv. For this reason, at first the shift of E a \ a dominates, but as rs increases the change in 
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E a \° M becomes more important, and thus crosses zero a second time for rs ~ 25 Sv, and is positive afterwards. 
This second zero will move to higher values of r$ as n increases, marking the position where the effects of rs and k 
compensate each other (not shown). The values of the shift are always relatively small in magnitude (< 0.05 Sv), as 
long as rs is not unrcalistically large. 

The dependency of on k indicates that an increased vertical diffusion does not affect M ov , but shifts L2 towards 
more negative values of E a , by stabilising the "OFF" state of the MOC and widening the ME region. This translates 
into the shift of towards negative values as k increases, as seen in Fig. [10] This stabilising effect of vertical mixing 
within the Atlantic basin for the "OFF" state was already recognised by Sijp and England [33]. In their study, they 
could not reach any permanently reversed state below a critical value of vertical diffusivity within the Atlantic basin. 
This is not the case in our study. This discrepancy is likely due to the use of perturbations too small to push the 
system to the "OFF" state with low diffusivity values in the work by Sijp and England [33J. Considering the values 
of vertical diffusion used in most numerical models (O(10~ 4 m 2 s -1 )), we may expect the skill of M ov in identifying 
the ME regime to be often low. The fact that M ov has been instead successfully used to identify the ME regime in 
different numerical models (e. g. |I1|6][32]) is a sign that a compensation between the effects of vertical diffusion and 
horizontal advection (and/or diffusion) is taking place in numerical models. 



When rj is increased, goes from positive to negative values, if rs, k or both are greater than zero (Fig. 10 
bottom). This is due to the proportionality between M ov and the MOC strength: increasing the hydraulic constant 77 
will increase M ov , requiring a higher E a for bringing it to zero, while the L2 position is independent of 77 (Eq. ([8])). A^ 
is always smaller than 0.1 Sv if rj is taken in a range giving MOC rates not too far from reality (77 ~ 2 — 6 x lO^ns" 1 ). 
If the gyre in the south and vertical diffusion are both set to zero, A^ is identically zero independently of rj, as 
discussed above. 



IV. SUMMARY AND CONCLUSIONS 



A box model for the MOC was developed, focusing on the effect of the exchanges between the Atlantic Ocean and 
Southern Oceans on the stability of the Atlantic MOC. The model includes a shallow box of variable depth, providing 
a basic representation of the Atlantic pycnocline depth dynamics in the tropics and subtropics, similarly to what 
was done in [14] and J07. The pycnocline depth is set by a balance between inflow from the Southern Ocean due to 
Ekman pumping, outflow due to baroclinic activity near the southern subpolar front, upwelling at low latitudes due 
to vertical diffusion and downwelling at the high northern latitudes. Differently from previous studies, the scaling 
for the downwelling flux in the north depends on the density difference between the northern North Atlantic and the 
region above the pycnocline north of the ACC, which is represented by a separate box. Furthermore, the transport 
(of salinity) between the ACC and the Atlantic is treated in a different way than in earlier studies: Ekman inflow and 
eddy outflow are considered separately for what concerns their associated salt transport. 

Similarly to what obtained in J07 and Gnanadesikan (14j . this configuration produces an inter-hemispheric MOC 
only if sufficiently strong wind stress is acting over the Southern Ocean. The MOC is, however, still buoyancy 
controlled, as changes in the north-south density difference, determined in the model only by changes in freshwater 
forcing or transport, can collapse the MOC. The main advance with respect to previous works, is the ability to 
reproduce the sensitivity of the MOC to changes in the freshwater transport by the southern subtropical gyre in the 
Atlantic Ocean. These results point to the fundamental importance of the region south of the tip of Africa and north 
of the ACC in determining the MOC stability. 

In particular, the freshwater transport at the latitude of the southern tip of Africa, by either the overturning or 
the azonal circulation, can change completely the response of the MOC to perturbations in the surface freshwater 
flux. In this view, the net freshwater import into the Atlantic basin by the meridional overturning circulation, M ov , 
is playing a fundamental role, being associated with the growth of perturbations due to the salt-advection feedback. 
If the salt-advection feedback is the dominant feedback connected with a MOC collapse, the sign of M ov completely 
determines whether a permanent collapse is possible. If other responses are important (gyre circulation, vertical 
diffusivity, atmospheric feedbacks or others not considered in this paper), the sign of M ov will not be a perfect 
indicator of multiple steady states anymore, its skill being dependent on the strength of the salt-advection feedback 
in comparison with other responses. For a parameter set representative of the real ocean, we can expect M ov to be 
well below 0.1 Sv when entering the ME regime. Considering the effect of noise and internal variability, this difference 
from zero is unlikely to be important. On the other hand, if vertical diffusion is providing a strong feedback during a 
MOC collapse (in particular in numerical ocean models with high vertical diffusivity — both physical and numerical) 
the skill of M ov as indicator of ME may be reduced. The observation that M ov actually is a good indicator in different 
numerical models [4, 6, 32 , suggests that a compensating mechanism between different feedbacks is operating, in 
particular between the response due to vertical diffusivity and horizontal advection by the southern subtropical gyre. 
Apart from its sign, a downward trend in M ov in response to a slowly varying external forcing |38j is a robust signal 
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of an approach to the collapse point of the MOC, independently of the source of the freshwater perturbation (i. e. any 
term of the freshwater budget Eq. Q). 

The box model enables to speculate on the importance of the latitude at which M ov is computed: "Why thirty 
degrees south?" or, in other words, why should M ov be computed at the southern edge of the Atlantic Ocean? One 
hypothesis is that the latitude of the southern tip of Africa marks the point north of which the MOC can be considered 
as a coherent flow with approximately constant water properties in numerical models. Southward of this latitude and 
down to the northern end of the ACC, the interaction between the Ekman inflow from the Southern Ocean and the 
salt transport from the north by the subtropical gyre determines the density of the waters that flow into the Atlantic 
basin. North of this point, salinity anomalies entering the Atlantic Ocean will eventually be advected to the northern 
downwelling regions where they will translate to a MOC strength anomaly (neglecting feedbacks other than the salt- 
advection feedback). Whether these anomalies are subject to a positive feedback loop or a negative one depends on 
the sign of the freshwater transport at the southern entrance of the basin. A second hypothesis, possibly related 
to the first one, concerns the presence of zonal boundaries north of the southern tip of Africa, which can support a 
zonal pressure gradient across the basin, controlling the MOC strength through planetary geostrophic balance [34] , 
Changes in this region may affect the entire basin, setting the boundary conditions for the MOC within the basin. 

We have demonstrated that the box model is a useful diagnostic tool to understand results from GCMs, and 
that M ov can be a good indicator of the stability of the MOC, but one may ask if these results are relevant to 
interpret the results from high-resolution GCMs and, more importantly, observations from the real ocean. In order 
for our conclusions to hold in those cases as well, the main assumption that must be satisfied is that the meridional 
overturning circulation should have a meridionally coherent response to freshwater perturbations on long time scales, 
at least within the Atlantic basin. Finally, even if the role of meridional density gradients in driving the MOC is 
corroborated by various studies (see e. g. [25, 35 ), the relevance of the scaling for the downwelling flux in an eddying 
ocean with complex geometry is still to be demonstrated. While the downwelling in the subpolar gyre seems to be 
well described by the theory of Spall |36j . the connection between the downwelling in the subpolar gyre and the 
overturning circulation and its scaling law on the global scale is still unclear. 
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Appendix A: Derivation of E a | m: ov and Ae 



To compute the E a value needed to bring M ov to zero with q$ > 0, the system ( 10 ) has to be solved. 
In the non diffusive limit k — > Eq. (lOe) reduces to qs — qN = 0, a second order equation in D, which can be 
solved giving as positive solution: 
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This result can be substituted into Eqs. (10a)-(10dl, and salinity can be eliminated giving: 

r N +r s + {T ts - T n ) (D\° M J 2 a V 
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= E a 



The latter is solved for E a giving E a \° M ^ in the case of n — s- 0: 
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(Al) 



(A2) 



(A2) can be combined with Eq. ([9| to give in the non diffusive limit, Eq. (11). 

The method outlined for the non diffusive limit can be used for the model including vertical diffusion as well. Also 
in this case, the mathematics involved is very simple on conceptual grounds, but the large expressions obtained render 
the problem tedious with pencil and paper. The final result is obtained with Mathematica software. Eq. (lOd), in 
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this case a third order algebraic equation, can be solved for pycnocline depth, and the positive solution reads: 
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After eliminating salinity, Eqs. (10a)-(10d| reduce again to (All, and E a \° M is still given by Eq. (A2), but with 
D\° Mn , and 

9JvIm „ f° r the fi n he vertical diffusion case. 
The difference between E a \ qN and E a \° M then gives A^, which can be written as: 
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In the limit of no vertical diffusion, q e \q N = qEk if qN = 0, and thus reduces to: 
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which can be written in the form of Eq. (11 1 when <7jv|m i s written explicitly. 
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Vtot 

A 

(3 
L x s 
fs 
Aqm 
tn 
V 



3.0 x 10 17 m 3 
1.0 x 10 14 m 2 
8 x lO^psu" 1 
3.0 x 10 7 m 

1700 m 2 s" 1 
5 x 10 6 m 3 s- 1 
3.0 x 10 4 ms- 1 



Po 

L X A 
T n 

E s 
rs 

K 



9 x 10 15 m 3 
1027.5 kg rrr 3 
5°C 
1.0 x 10 7 m 
5.0°C 
0.25 x 10 6 m 3 s- 1 
10 x 10 6 m 3 s- 1 
m 2 s~ 



Vn 

a 
So 

Li> 

T ts 
E„ 

T 



3 x 10 15 m 3 
2 x lO^K- 1 

35 psu 
1.0 x 10 6 m 
10.0°C 
-0.1 x 10 6 m 3 s- 1 
O.lONm- 2 



Table I. Reference values of the parameters used in Eqs. |2| and |3|. 
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696 m 


S n 


35.0 psu 


s d 


s t 


35.2 psu 


Sts 


34.6 psu 


S s 




18.8 Sv 


qs 


17.4 Sv 





35.0 psu 
34.4 psu 
1.4 Sv 



Table II. Steady solution for the reference configuration of the model (using the parameters in tab. |D. 




Figure 1. Sketch of a bifurcation diagram for the MOC strength as a function of a generic freshwater forcing strength 7. The 
blue lines mark the two stable steady state solutions, "ON" and "OFF"; the grey line marks the unstable solution connecting 
the two stable ones. At points L\ and L2, marked by green circles, the stability of the steady state solution changes. The 
range of 7-values between the two dashed vertical red lines is the ME regime, where two stable states coexist under the same 
boundary conditions. Changing 7 further, the system jumps from one solution to the other at the two limit points L\ and L2, 
as indicated by the red arrows. 
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Figure 2. Structure of the box model, with flow pathways connecting the different boxes. Red arrows represent net volume 
fluxes (with names in the same colour), green arrows represent gyre exchanges between the different boxes (with names in the 
same colour), blue arrows are water vapour transports through the atmosphere (with names in the same colour). In black, the 
boxes names. 
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Figure 3. Sensitivity of D (top panel, blue), qs (lower panel, blue) and <?jv (lower panel, red) to the wind stress r. All other 
parameters are kept at the reference values of Tab. [I] 
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Figure 4. Sensitivity of D (top panel, blue), qs (central panel, blue), gjv (central panel, red), M ov (lower panel, blue) and M az 
(lower panel, red) to the gyre exchange between the boxes ts and t (r\s). All other parameters are kept at the reference values 
of Tab. H 
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Figure 5. Sensitivity of D (top panel, blue), qs (central panel, dashed) and gjv (central panel, full) to the symmetric freshwater 
flux E s . The different colors (red, green and blue) refer to different r\ values: (1.5,3.0,4.5) x lff'ms^ 1 respectively. 
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Figure 6. Sensitivity of D (top panel, blue), qs (central panel, dashed), gjv (central panel, full) and M ov (lower panel) to the 
asymmetric freshwater flux E a . The different colors (red, green and blue) refer to different 77 values: (1.5,3.0,4.0) x 10 4 ms _1 
respectively. L\ marks the position of the saddle node bifurcation ending the "ON" state of the MOC. After Li, the solution is 
unstable, and for higher values of E a only a collapsed state of the MOC is possible. L2 marks the position of the saddle node 
bifurcation ending the "OFF" state of the MOC. The plot of M ov stops before the saddle node bifurcations for clarity. 
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Figure 7. Continuation of the limit point Li, determining the collapse of the MOC, (see Fig. [6|. A stable "ON" state exists 
only below the lines, marking the position of L\. On the top panel, the continuation in (E a ,rs) is shown, while in the lower 
panel the same regime diagram is shown using M az instead of rs for plotting. The different colors (red, green and blue) refer 
to different r\ values: (1.5,3.0,4.5) x 10 4 ms -1 respectively. 
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Figure 8. Changes in the bifurcation diagram of the MOC as a function of E a , as the strength of the gyres are changed; the full 
(dashed) line refers to qN {qs)- On the top panel, rs is changed, while in the lower panel tm is changed. The different colors 
(red, green and blue) refer to different r$ or rjv values: 5, 10, 15 Sv respectively. For the case of r^, the discontinuity in the 
MOC strength at the point where qs = is due to the fact the when qs reverses the sensitivities of the MOC change abruptly. 
Li and L2 mark the limit points ending the "ON" and "OFF" state respectively. 
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Figure 9. Continuation of the limit point L\ determining the collapse of the MOC (see Fig. [fjj. A stable "ON" state exists 
only below the lines, marking the position of L%. On the top panel, the continuation in (E a ,rN) is shown, while in the lower 
panel the same regime diagram is shown considering q$ instead of E a , that is the strength of the flux in the south at the point 
of collapse. The different colors (red, green and blue) refer to different r\ values: (1.5,3.0,4.5) x 10 4 ms _1 respectively. 
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Figure 10. Dependence of Ae on rs (top), k (centre) and r\ (bottom). A_e is shown on the top (central) panel as a function rs 
(k), keeping k = Om 2 s _1 (rs = OSv); all other parameters are kept at the reference value of Tab.|l] In the bottom panel, A_b 
is shown as a function of rj with all other parameters as in Tab. [I] 



